Volovik effect in the ±s-wave state for the iron-based superconductors 
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We studied the field dependencies of specific heat coefficient y(H) = linij ^qC(T \H) /T and thermal con- 
ductivity coefficient limJ'_ > o K (^,//)/7 , of the ±s-wave state in the mixed state. We found that it is a generic 
feature of the two band s-wave state with the unequal sizes of gaps, small A$ and large A^, that Doppler shift 
of the quasiparticle excitations (Volovik effect) creates a finite density of states, on the extended states outside 
of vortex cores, proportional to H in contrast to the \f~H dependence of the d-wave state. Impurity scattering 
effect on the ±s-wave state, however, makes this generic //-linear dependence sublinear approaching to the 
\fH behavior. Our calculations of lim7-_>o k(T,H)/T successfully fit the experimental data of Ba(Fei_ v Co x )2 
As2 with different Co-doping x by systematically varying the gap size ratio R = |As|/|Al|. We also resolve 
the dilemma of a substantial value of y(H — > 0) but almost zero value of UmT^oK(T,H — > 0)/T, observed in 
experiments. 
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Introduction: As to the pairing symmetry of the recently 
discovered Fe-pnictide superconductors (SC) id, the its- 
wave state 12ll3[ is considered as the most natural and promis- 
ing pairing state. Genuinely being an s-wave gap state, this 
model is expected to provide exponentially activated behav- 
iors of various thermodynamic and transport properties due 
to full gap(s) around Fermi surface(s) (FS) in the supercon- 
ducting phase. Angle resolved photoemission spectroscopy 
(ARPES) experiments [4] and early measurements of penetra- 
tion depth of M-l 1 1 1 (M=Pr, Nd, Sm) fl and (Ba,K)Fe 2 As 2 
|0] confirmed this expectation. However, some power laws of 
temperature dependence measured, for example, with the nu- 
clear spin-lattice relaxation rate T^ 1 ~ I 3 |J as well as the 
penetration depth X(T) ~ f 2 -- 2 - 5 on Ba(Fe,Co)2As2 10] were 
not consistent with this simple picture. But these power law 
dependencies were successfully explained theoretically with 
the ±s-wave pairing state by the sign-changing feature of the 
order parameter (OP) and the related noble impurity effect 
and hence strengthened the status of the ±s-wave state 
as the pairing symmetry of Fe-pnictide superconductors. 

However, recent measurements of specific heat coefficient 
y(H) = C(H) /T and thermal conductivity coefficient k(H) /T 
in the mixed state and their field dependencies are posing a 
new challenge to the ±s-wave pairing model. Several mea- 
surements of specific heat coefficient showed a strong field 
dependence: Ba .6K .4Fe2As 2 UM ( x H), (Feo. 92C oo.oR)? As?. 
I0 (sublinear in H), and LaO .9F i_ 5 FeAs (°c y/H). 
And thermal conductivity measured on Co-doped Ba-122 111411 
and K-doped Ba-122| 15] also showed various field dependen- 
cies. In particular, Tanatar et al. 11411 measured k(H)/T of 
Ba(Fei_ x Co x )2As2 for several values of x and found that the 
field dependence of k(H) /T continuously evolves from expo- 
nentially flat, near linear in H, and to near y/H behavior with 
increasing Co doping x. 

Another puzzling observation is that, while the zero field 
limit of y(H = 0)/y nurma i in several experiments II 1 ll — 1 1 311 
showed a substantial values, the zero field limit of k(H = 
0)/T ffl [l5tl of the same compounds with similar dopings 



approaches to negligibly small value. It implies that there ex- 
ist a large amount of zero energy excitations but they do not 
contribute to thermal conductivity. Neither a nodal gap state 
nor a simple isotropic s-wave gap state can be easily compat- 
ible with these observations. 

In this paper, we study the field dependence of the spe- 
cific heat and thermal conductivity coefficients of the ±s-wave 
model in the mixed state using a semiclassical approximation 
which causes Doppler shift of quasiparticle excitations outside 
of vortex cores (so-called "Volovik effect" Recently, 
Mishra et al. 11711 studied the field dependence of thermal 
conductivity on the isotropic ±s-wave model with equal size 
gaps using a different method [18] and concluded that this 
model is incompatible with the experiments HUH]- Their 
result agrees with ours in that particular case. In the current 
paper, however, we show that it is essential to take into ac- 
count the unequal sizes of gaps in the ±s-wave model. This 
feature is not only more realistic but also introduces genuinely 
new physics to explain experimental data. We found that the 
unequal size of two isotropic s-wave gaps causes the field de- 
pendencies of k(H)/T as well as y(H) to be linear in H at the 
low field regime. 

This generic //-linear dependence in the ±s-wave model is 
in contrast to the \[H dependence in the d-wave state iflfill al- 
though both are due to Doppler shift effect on the extended 
states outside of vortex cores - it should also be contrasted to 
the //-linear contribution to specific heat from the core bound 
states in a single s-wave state [19] that cannot contribute to 
thermal conductivity. We then show that the impurity scatter- 
ing, combined with Volovik effect, make the generic //-linear 
dependence gradually modified to become sublinear and to 
approach to the \Jli dependence. 

Formalism: For the generic two band model of the ±s-wave 
pairing state, we assume two isotropic s-wave order parame- 
ters (OPs) As and Al of opposite signs on the two representa- 
tive bands of the Fe pnictide materials. As and Al represent 
a smaller and a larger size gap on the corresponding bands, 
respectively. All energy units in this paper is normalized by 



2 



At. 

We use the semiclassical approximation where the vector 
field term in Hamiltonian | A(r) • k is replaced by the circulat- 
ing superflow velocity as v s (r) • k that causes Doppler shifting 
of the quasiparticle excitations. The validity of this method 
is discussed in length in Ref. lll7ll and it is a reliable approxi- 
mation from the low to intermediate field regime but could be 
questionable approaching H c i, Our results, therefore, should 
be viewed with reservation at the high field regime where a 
better method 111 811 should be applied in principle. 

The single-particle Green's function of band a(= S,L) in 
Nambu matrix form including Doppler shift of the quasiparti- 
cle excitations is a function of r the distance from the vortex 
core in addition to the usual momentum k and frequency CO as 
follows |d. 



G a (k,r,co) = 



[CO + v s (r ) ■ k] Tq + ^ a (k) T 3 + A a X 1 

[co + v s (r).kP-^ (fe) _ A 2 



(1) 



where x, are Pauli matrices and £, a (k) is the quasiparticle en- 
ergy of band a. From the above Green's function we obtain 
the local density of states (DOS) of each band asN a ((d,H,r) = 
— ^Tiim^ l(; G a (k,r,co). The Doppler shifting energy is given 
as v s (r) • k = — - cos = b-£ cos 8 with normalized distance 
p = rl\ (t, = coherence length) and "b" a constant of order 
unity. 

We use the above Green's function for 1 < p < Rh/^ with 

magnetic length Rh = < ^\f^ i'&Q a f ux quanta, H magnetic 
field, and "a" geometric factor of order unity) because the 
core region (p < 1) is thermodynamically negligible if fields 
is not too close to Eq.(l) shows that the gaps collapse 
whenever the Doppler shifting energy v s (r) k becomes larger 
than the gap energy A„. The larger gap At does not collapse 
even with the largest Doppler shifting of v s (p = 1) ■ kp ~ At 
at the boundary of vortex core. But the small gap As col- 
lapses for the region of 1 < p < to create a finite DOS 

N s ((0 = 0,H,r) « jV™"' M/ . With this observation, the mag- 
netic unit cell averaged DOS N a ((0,H) =<N a ((d,H,r) > C ell= 
f£ H dr 2 N a ((i),H,r)/KRj 1 is readily obtained at co = as fol- 
lows. 
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Above Eq.(3) holds as far as As < At and shows that 
Volovik effect immediately create a finite DOS with the 
isotropic ±s-wave state and there is no threshold value of 
magnetic field H* to collapse the small gap A$. Its generic 
field dependence is linear in H and its slope is proportional 
to ps (s^) 2 - We will show later with numerical calculations 
that impurity scattering will smoothen this generic linear-in- 
H field dependence and make it more sublinear and closer to 



<*■ \/H. Therefore the impurity effect is important to under- 
stand experiments. 

The impurity scattering is included by the 1 -matrix 
method, suitably generalized for the ±s-wave pairing model 
ifioll . And impurity induced selfenergies renormalize the fre- 
quencies and OPs as co -» co = co+E^co) + E° (co), and A s .l 
h,L = A sx + E s (co) +E|(co), with Esi(co) 



where T = -^f-', ni mp the impurity concentration and N to t = 

n s+ N l is the total DOS - The T -matrices T ' 1 are the Pauli 
matrices x ' 1 components in the Nambu space. We refer read- 
ers to Ref.[10] for more details and it is straightforward to 
incorporate the local Green's function Eq.(l) with Doppler 
shifting into the T -matrix method. 

After calculating the averaged N a ((o,H) for all frequencies, 
specific heat is calculated as 
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Similarly, thermal conductivity is calculated with i20ll 
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where z = CO + v s (r) • kp. And then longitudi- 
nal and transversal thermal conductivities are cal- 
culated as K||(r,H) = J cell d 2 nc(T,H,r)/%R 2 H and 
K~l l (T,H) = J cell d 2 rK- 1 (T,H,r)/'KRj I , respectively. 

Numerical results and discussions: Since we consider two 
band model, total specific heat and thermal conductivity are 
the sum of two contributions from each band. In order to de- 
termine the precise contributions, we need to know the normal 
state DOSs N$ L as well as Fermi velocity Vp(s,L) °f eacn band. 
However, in all numerical calculations in this paper, the large 
gap band is almost gapped even with a finite amount of impu- 
rity concentration so that its contribution to specific heat and 
thermal conductivity are very small compared to the contribu- 
tion from the small gap band. Therefore it is not very mean- 
ingful to determine the precise weighting of each band and 
we simply put equal weighting to each band for convenient 
illustrations. For the field dependence of the OPs, we use a 
phenomenological formula A$.l(H) = A a SL yJ\ —H/H C 2- 

In Fig.l, we show numerical calculations of a pure case 
(r/A° = 0.001) with a typical gap size ratio |A°/A°| = 0.3. 
Figures 1(a) and (b) show the magnetic unit cell averaged 
DOSs Ns.l(u>,H) of the small and large gap bands, respec- 
tively, as a function of magnetic field H. As expected, the 
small gap DOS Ns((0,H) immediately starts collapsing with 
field and reaches a constant value Nf udl at CO = for H > 
0.\H c2 (see Fig. 1(c)) while the larger gap DOS N L ((0,H) 
remains gapped until H — > H C 2- However, the gapped re- 
gion in Nl((0,H) in frequency axis become extremely nar- 
row with increasing field which will affect specific heat coeffi- 
cient at high fields (see Fig. 1(d))). Fig. 1(c) shows normalized 
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FIG. 1: (Color online) (a-b) Normalized DOSs with Doppler shifting 
by magnetic field H, A^(co,/f) and Ni(a>,H) for small gap and large 
gap bands, respectively, in clean limit (r/A^ = 0.001). (c) Normal- 
ized DOSs at co = 0, Ns,L,tot(® = 0.H)/N' t , ° t r " w ' of small gap and 
large gap bands, and the total, respectively, (d) Normalized specific 
heat coefficients Js.l.ioi {H) Irfot" 1 " 1 °f sma 'l g a P ar >d large gap bands, 
and the total contribution, respectively. 



zero energy DOSs Ns,l(u> = 0,H) and their total N tot (0,//) = 
Ns{(ti = 0,H) + A\(co = 0,H). As explained above, beyond 
H s» 0.1H C 2, the small gap A^CO = Q,H) reaches its full value 
(= 0.5N® ul ) and the large gap iVz,(co = 0,H) remains zero until 
very close to H C 2- Fig. 1(d) shows normalized specific heat co- 
efficients Ys,L,tot(H), respectively. Some difference between 
Nl(Q,H) and Jl(H) at the high field region is due to the finite 
temperature (we used T /A£ = 1 /50 to calculate the T — » 
limit in this paper) and very narrowly gapped Nl( CO,//) at the 
high field region as seen in Fig. 1(b). 

In Fig. 2(a) and (b), we show calculations of total thermal 
conductivity and specific heat coefficients with varying gap 
size ratio \A° S /A° L \ = 0.2,0.3,0.5,0.7, and 0.9, respectively. 
Here we used impurity scattering with F/A^ = 0.05 and uni- 
tary impurity (c = 0) both for intra- and interband scatter- 
ing. These values are not particular choices for the plots. A 
wide range of impurity concentrations F and different strength 




FIG. 2: (Color online) (a) Normalized total transverse thermal con- 
ductivity coefficient lim^o [ K x,ror(//) /r]/[Kjv., 0/ /T] for different 
gap size ratios, |Ag/A£| = 0.2,0.3,0.5,0.7, and 0.9, respectively, (b) 
Normalized specific heat coefficient y tot (H) /Y/o/jv f° r different gap 
size ratios as in (a). All calculations include the same concentration 
of impurities T/Ai = 0.05 with unitary scattering limit (c = 0) both 
for intra- and interband scattering. 

(c 7^ 0) of impurity scattering produce qualitatively similar be- 
haviors. Fig. 2(a) shows normalized total transverse thermal 
conductivity lim T ^ Q [K ± , s (H)/T + K^{H)/T}/[y^ 1 /T]. 
Only the transverse K±/T (where the external field is applied 
as H || c and the thermal current is measured in the plane as 
j,/, 1 1 ab) are shown because those are the measured thermal cur- 
rent in most of experimental settings ffl fl5ll . The longitudi- 
nal thermal conductivity coefficient Ky /T (where j,/, || c) be- 
haves a bit more concave, in particular at the low field regime, 
because it corresponds to a parallel circuit of resistors while 
the transverse one corresponds to a series circuit. At the 
higher field regime they become indistinguishably similar to 
each other. 

The main feature of the results in Fig. 2(a) is the systematic 
evolution of the slope of K±(H) /T at the low field limit. For 
a large gap size ratio of |A§/A£| = 0.9, lim//^o Kx [H)/T be- 
comes exponentially flat similar to the behavior of a single gap 
isotropic s-wave state. Then with decreasing gap size ratio of 
|A§/A£|, the low field slope quickly increases and the overall 
H dependence of Kx (H) /T becomes concave down and close 
to the v 7 // behavior for |A§/A£| = 0.2. This concavity can be 
made even stronger with decreasing the gap size ratio | A^/A£| 
and increasing impurity concentration F. 

Another important feature is that the values of K±(H)/T 
(also K|| (H)/T although not shown here) in the zero field limit 
are negligibly small for all cases despite substantial impurity 
scattering; this is even true with the |A^/A£| = 0.2 case which 
shows the behavior ps \/H as in the d-wave case. In fact, these 
extremely small values of the thermal conductivity coefficient 
K±(H)/T in the zero field limit was argued as an evidence 
of an isotropic s-wave gap nature II 1 411 . However, it is a very 
puzzling feature when we note that the several experiments 
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1 11-3 observed substantial values of the specific heat co- 
efficients y(H —> 0) with the same compounds with similar 
dopings. Surprisingly, these seemingly conflicting behaviors 
are also obtained in our theoretical calculations. In Fig. 2(b), 
the normalized specific heat coefficient y(H — > 0) /jtot,N, with 
the same parameters as in Fig. 2(a), show substantial values in 
the zero field limit (0.5y ror jv is the full value of the small gap 
band contribution). 

It is a common knowledge that the DOS near Fermi level 
similarly contributes to both specific heat and thermal con- 
ductivity. However, the discrepancy between y(H —> 0) and 
k(H — >• 0)/T is possible due to the coherence factor of super- 
conductivity. The kernel of thermal conductivity (see Eq.(6)), 
being an energy current-energy current correlation function, 
has a destructive coherence factor ("— " sign in the numera- 
tor of the last term in Eq.(6)) and the specific heat does not 
have such coherence factor. In the case of d-wave pairing, the 
same coherence factor becomes very weak in the low energy 
limit because the nodal gap Ad(0) linearly disappears, so that 
y(H — s- 0) and k(H — » 0)/T behave similarly. 

As to the field dependence of j(H), with a relatively large 
gap size ratio, such as | / A£ | = 0.7 and 0.9, y(H) is very 
linear in H for a substantial region of fields (up to « H C 2/2) 
despite a finite y(H —> 0). Decreasing the gap size ratio to 
\A° S /A° L \ = 0.5, 0.3 and 0.2, the field dependence of y(H) be- 
comes gradually more concave down. This behavior is in ex- 
cellent agreement with the measurements of Bao gKo.4Fe2As2 
iflHl (« H), (Feo.92Co .o8)2As 2 10 (sublinear in H), and 
LaO .9F .i- S FeAs Q (« 

Conclusion: Using semiclassical approximation, we have 
calculated specific heat coefficient y(H) and thermal con- 
ductivity coefficient k(H)/T of the ±s-wave model in the 
mixed state with including impurity scattering. We found that 
Doppler shift on quasiparticle excitations on extended states 
immediately induce zero energy excitations proportional to H 
from a smaller gap band. There are no threshold values of 
small gap size A* s and magnetic field H* to induce such zero 
energy excitations. This steep increase ofy(H) and k(H)/T 
linear in H is a generic behavior of any two gap isotropic 
s-wave superconductors with unequal sizes of gaps. There- 
fore, our results also provide natural explanation for the ther- 
mal conductivity measurements in MgB2 [21] and NbSe9 H22ll . 
The sign-changing feature of the ± s-wave state is not a pri- 
mary source of the field dependence of y(H) and k(H)/T 
but plays an important role through impurity scattering. The 
impurity scattering effect with sign-changing OPs II 1 CJfl ef- 
ficiently accumulate zero energy impurity band near Fermi 
level and reflects it to y(H —> 0) (see Fig. 2(b)) and makes the 
field dependence of k(H)/T and y(H) more concave down 
and smoothly curved. 

Having found substantial contribution of low energy exci- 
tations from extended state outside of the vortex cores, we 
justified the ignoring of the core region contribution a pos- 
teriori. We compare the numerical calculations to the ther- 
mal conductivity measurement [14] of Ba(Fei_ T Co T )2As2 and 
showed that the systematic evolution of the thermal conduc- 



tivity coefficient k(H)/T with Co-doping x is well explained 
by changing the gap size ratio | A° /A£ | in our theoretical cal- 
culations, which is also consistent with the Fermi surface evo- 
lution in Ba(Fei_ r Co^)2As2 with electron doping by Co in the 
rigid band picture. We also showed that the conflicting behav- 
iors between j(H) and k(H)/T in the H — > limit is a nat- 
ural consequence of the superconducting coherence factor of 
thermal conductivity kernel. In summary, the ±s-wave state 
with isotropic s-wave gaps of unequal sizes can consistently 
explain the field dependence ofy{H) llTTl - [l3ll and k(H)/T ob- 
served in Co-doped Ba-122 Hj and K-doped Ba-122 111. 
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